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ABSTRACT 


Standard  reliability  prediction  formulas  lot  multi-component  .y.Un.  n*o 
the  assumption  of  statistical  Independence  of  the  conditions  of  the  conponents 
after  a  fined  period  of  tine  in  a  hypothesised  environment.  Althouph  labora¬ 
tory  and  field  experience  shoos  that  this  assumption  is  not  stays  valid  It 
persists  ns  a  basis  for  reliability  modeling  among  practitioners  of  the  art. 

The  most  relevant  reason  for  this  Is  that  engineers  -morally  arc  not  noted 
for  their  knowledge  of  the  mathematics  of  probability  and  one  soon  discovers 
that  the  assumption  of  statistical  independence  among  components  usually  leads 
to  the  Simplest  mathematics.  A  second  reason  why  the  Independence  assumption 
is  so  predominant  is ‘that  engineers  and  analysts  are  not  clear  on  what  alterna¬ 
tives  should  be  pursued  or  oven  what  the  alternatives  night  be  so  there  uould 
be  little  point  In  maklnp,  assumptions  the  modeling  implications  of  which  ate 
simply  not  understood.  There  is  .another  Justification  for  the  statistical 
independence  assumption,  that  being  the  fact  that  such  an  assumption  yields 
models  that  can  provide  bounds  on  system  reliabilities.  Analytical  models 
which  attempt  to  account  for  environmental  effects  on  component  failure  rates 
show  that  the  statistical  independence  assumption  often  leads  to  gross  over¬ 
estimates  or  underestimates  of  system  reliability.  Evans  [1]  argues  that  it 
is  almost  never  the  case  that  this  assumption  is  correct,  but  also  points  out 
that  other  models  which  attempt  to  quantify  or  correct  this  error  have  their 
own  limitations.  If.  on  the  other  hand,  it  could  bo  shown  that  an  alternative 
model  even  though  an  approximation  Is  indeed  a  closer  representation  to  the 
true  state  of  the  world  and  is  at  the  same  time  practical  to  use  In  a  computa¬ 
tional  sense  then  it  should  be  used.  The  question  of  validity  lies  with  the 
basic  assumptions  underlying  the  model  and  not  with  the  mathematics  Itself. 
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If  a  set  of  assumptions  alternatlva  to  th«  assumption  of  statistical 
independence  are  accepted  as  more  closely  characterizing  an  environmental 
situation  then  the  resulting  mathematical  predictions  would  be  tentatively 
preferred,  subject  to  experimental  verification. 

This  paper  reviews  some  previous  explorations  of  the  question  of 
statistical  dependence  of  operating  components  and  the  effects  on  system 
reliability  and  presents  some  results  not  previously  worked  out*  In  particu¬ 
lar  there  is  shown  to  exist  a  connection  among  three  models  that  do  not  re¬ 
quire  the  assumption  of  statistical  Independence  of  components.  Under  certain 
conditions  these  models  can  generate  the  same  reliability  prediction  even 
though  they  appear  to  be  derived  under  different  assumptions. 
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CHAPTER  I 


DEFINITIONS 

For  what  follows  it  is  necessary  to  distinguish  among  the  notions  of  the 
physical  dependence  of  two  events  (or  two  random  variables) ,  the  conditional 
stochastic  independence  of  two  events  given  a  third  event,  and  the  unconditional 
stochastic  independence  of  two  events. 

Definition  1  If  two  events  A  and  B  are  unconditionally  statistically  in¬ 
dependent  then  P(AB)  ■  P(A)  P(B), 

The  question  of  the  conditional  statistical  independence  of  two  events 
A  an'd  B  given  a  third  event  E  arises  when  components  arc  placed  into  a  random 
operating  environment.  The  events  A  and  B  represent  the  conditions  of  two 
components  after  a  suitable  period  of  tine  when  operated  in  the  presence  of  a 
randomly  selected  environmental  condition  E. 

Definition  2  If  two  events  A  and  B  jointly  conditioned  with  respect  to  an 
event  E  are  statistically  independent  under  F.  then 

P(AB/E)  *  P(A/E)  P(B/E). 

Let  ■  (E^s  i  e  1}  be  a  collection  of  mutually  exclusive  and  exhaustive 
events  called  an  environmental  profile.  If  A  and  B  arc  statistically  indepen¬ 
dent  under  the  environmental  profile  £  »  thcr*  P(A,B/E^).  **  P(A/E^)  P(B/Ej.)  for 
all  E^  e  £  .  Thus  if  AB  occurs  only  in  the  presence  of  an  E^  e  £  and  A  and 
B  are  statistically  independent  under  all  E i  e  £  ,  then" 

/  P(AB/i;.)dP(E.)  =  /  P(A/E  )  P(B/E.)dP(Ei). 

Ei  Ei 
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If,  in  addition,  A  and  »  are  unconditionally  statistically  independent 

then 

/  P(A/El)  P(B/E  )<1P(Ej)  -  /  P(A/Ei)dP(Ei)  /  P (3/lu±)dP (E±)  -  P(A)  P(B) 
Ki  Ei  ni 

The  intuitive  notion  of  the  physical  independence  of  two  components  is 
that  the  condition  of  each  component  neither  affects  nor  is  affected  by  the 
condition  of  the  other.  Following  Evans  (p.  343)  this  idea  is  restated 
quantitatively  as  follows. 

Definition  3  Two  events  A  and  B  are  physically  independent  under  an  environ 

mental  profile  £  if  and  only  if  they  are  statistically  independent  under  £  , 

i  p  that  ^ 

P(A,B/E1)  -  PU/E^  PtB/E^  for  all  Ej  c  £  . 

Thus,  if  the  malfunction  of  component  i  necessarily  destroys  component  j, 
intuition  states  that  they  are  certainly  not  physically  Independent  under 
such  an  environmental  circumstance.  If  such  an  event  were  Included  in  £  as 
Ej  with  positive  probability  mass  assignment  then  a  calculation  will  also 
show  that  the  events  A  and  B  representing  the  conditions  of  the  components 
are  not  statistically  Independent  under  E^  and  hence  not  physically  indepen¬ 
dent  under  E^. 

The  next  definition  characterizes  the  notion  of  "associated  random 
variables"  Introduced  by  Esary,  Proschan,  and  Walkup  [2]  in  1966.  This 
definition  is  relevant  to  what  follows  because  Pollyak  [3]  showed  in  1962 
that  under  a  random  environmental  profile  as  specified  by  Definition  2  the 
individual  component  reliability  functions  are  random  variables  which  under 
certain  conditions  are  associated  in  the  sense  of  Esary. 

Definition  4  Random  variables  X^,  X2,  ...  ,  ^  are  associated  if 
Cov  [U(X),  V(X)]  £  0  for  all  pairs  U,V  of  bounded  continuous  nondecreasing 
*It  is  also  the  case  that  if  P(AB/E1)  ■  P(A/Ei)P(B/Ei)  and  P(A/Et)  -  P(A)  for 
all  Et  then  P(AB)  -  P(A)P(B).  2 


functions  where  X  -  <Xl#  X2,  ...  ,  XR).  Esary,  at  si.,  showed  that  association 
has  the  following  proparties i 

Property  1  Any  subset  of  a  sat  of  associated  random  variables  is  also 
a  set  of  associated  random  variables. 

Property  2  If  two  sets  of  associated  random  variables  are  independent 
of  one  another  then  their  union  is  a  set  of  associated  random  variables. 

Property  3  The  set  consisting  of  a  single  random  variable  la 
associated. 

Property  4  If  X  ,  Xj,  ...  ,  X^  are  associated  then  any  set  of  non- 
decreasing  functions  S^X),  S2(X),  ...  ,  SM(x)  are  associated. 

•  The  concept  of  equally  ordered  functions  will  be  of  use  and  is  specified 
by  the  following  definition. 

Definition  5  Let  F(t,X)  and  G(t,X)  be  functions  of  a  parameter  t  and  a 
vector  argument  X  ■  (X^,  ...  ,  X^  where  Xt  is  a  real  number.  F  and  G  are 
said  to  be  equally  ordered  with  respect  to  the  variables  X^  ...  ,  X^  if 
for  any  values  of  X^  and  X2  of  the  vector  argument  X  the  inequality 

[FU.xp  -  F(t,X2)][G(t,X1)  -  GCt.Xj)]  >  0 

is  satisfied.  Monotonic  nonincreasing  or  nondecreasing  functions  of  a  single 
scalar  argument  defined  on  the  real  line  belong  to  this  class  as  well  as 
certain  convex  functions. 

If  F  and  G  in  Definition  5  are  functions  of  a  random  vector  X  -  (X,  ...  , 
Xj^),  then  the  variables  X1>  ...  ,  Xk  are  associated  in  the  sense  of  Definition  4. 
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CHAPTER  II 


THE  ENVIRONMENTAL  PROFILE  MODEL  OF  RELIABILITY 


It  is  evident  that  most  components  of  any  system  function  within  an 
environment  that  is  "uncertain"  in  the  sense  that  the  intensity  of  its 
defining  variables  exhibit  statistical  fluctuation  over  time.  Examples  of 
environmental  variables  are  temperature,  pressure,  humidity,  vibration 
frequency,  salt  content,  and  radioactive  emission  level. 

Let  the  environment  be  characterized  by  a  random  vector  X  ■  (X^,  ...  ,  X^ 
where  Xi  represents  the  i-th  environmental  variable  defined  over  a  suitable 
range  of  values.  The  environmental  profile  is  specified  by  the  set 

^  “  {(X^,  ...  ,  X^) !  a±  <  X±  <  8^,  i  “  1,  ...  ,  k)  . 

Let  the  joint  probability  density  of  X  be  u(X)  ■  u(X^,  ...  ,  X^) 

Let  R  (t,X)  be  the  conditioned  reliability  of  the  1-th  system  component 
where  t  denotes  tine  and  X  indicates  that  reliability  is  a  function  of  (or 
conditioned  upon)  X.  Pollyak's  model  of  component  reliability  is  obtained 
as  an  average  over  the  environmental  profile  £  I 

R(t)  -  /  R(t,X)  dw(X)  ,  Xef,  (1) 

X 


The  reliability  of  n  components  in  series,  assuming  mutual  statistical 
independence  of  the  components  under  £  is 


/  IT  R  (t,X)  du(X) 


X  i 


Xe£. 


(2) 


R^^X)  is  usually  assumed  to  be  monotonic  nonlncreas-ing  in  t  but  what  Is 
more  Interesting  is  the  effect  on  Fs,(t)  of  averaging  the  product  of  the  component 

reliabilities  over  the  sample  space  of  environmental  possibilities.  Pollyak 
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v  rod  this  question  by  provin;;  the  following  theorem. 


theorem  1  Let  i^Ct.X)  and  Hj(t,X)  be  functions  of  the  k-dimensional  vector 

a  rumen t  X  that  are  equally  ordered  in  the  sense  of  Definition  5.  Let  X 

belong  to  a  set  £  and  let  u>(X)  be  a  weight  function  defined  on  &  which  is 

everywhere  nonnegative  and  J  w(X)uX  ■  1. 

X 

Then 

f  R,  R.  u(X)dX  -  Jr.  *(X)dX  J  R,  u(X)dX  •  0  O) 

xU  1  d  » e  **  J 

Equation  (3)  implies  that  Cov(Ri,Rj)  i  0.  Consequently  the  effect  of  im¬ 
bedding  component  reliabilities  within  a  random  environment  and  averaging  over 
t;,e  environmental  possibilities  with  respect  to  a  distribution  to  obtain  the 

unconditional  reliability  is  the  following: 

(a)  For  serial  systems  reliability  is  at  least  as  high  as  the  value 
predicted  by  making  the  assumption  of  unconditional  statistical  indeper.den  _e 
\,i  h  respect  to  the  failures  of  the  components. 

(b)  For  parallel  systems  reliability  is  at  least  as  low  as  the  value 
predicted  by  making  the  assumption  of  unconditional  statistical  independence 

w.ong  the  conditions  of  the  components. 

In  other  words,  if  the  environmental  profile  model  as  defined  by  Equation 
(2)  is  a  more  realistic  reflection  of  physical  reality  and  if  component  failure 
laws  satisfy  conditions  of  Theorem  1  then  standard  prediction  formulas  for 
serial  reliability  give  underestimates  and  for  parallel  system  reliability  they 
ive  optimistic  results.  (This  interpretation  follows  since  pairs  of  products 
of  equally  ordered  functions  are  equally  ordered.  For  an  alternate  justifica¬ 
tion  see  Theorem  2.) 

It  seems  plausible  that  the  assumption  of  nonotonicity  of  the  component 
failure  law  with  respect  to  the  degree  of  severity  of  the  environment  X 
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•'i  uld  be  satisfied  in  most  problem  situations.  Tne  environmental  variables 

id  be  represented  either  directly  or  implicitly  in  terms  of  the  moments  of 

—  —  -0.1 1  — 
r..u:  iti(c .X) .  For  example  if  R^t.X)  -  e  then  0^  depends  upon  X  perhaps 

as  a  linear  regression  upon  (X^,  ...  ,  X^)*  Equation  (3)  does  not  imply  that 
on a  needs  to  imbed  components  in  a  random  operating  environment  which  pro¬ 
duces  a  positive  correlation  among  component  conditions  in  order  to  increase 
system  reliability.  To  attempt  to  do  so  to  an  extent  greater  than  what  is 
usually  the  case  anyway  may  impair  system  reliability.  Equation  (3)  does 
imply  however  that  components  operating  in  a  random  environment  and  having 
failure  lavs  that  are  monotonlcally  related  to  the  environmental  parameters 
t  :nd  to  fail  together  or  function  properly  together.  That  is,  if  the 
environmental  model  applies,  knowledge  that  a  component  has  either  failed 
(cr  survived)  increases  the  betting  odds  that  another  component  also  failed 
survived). 

The  fact  that  the  environmental  profile  model  is  a  more  accurate  re¬ 
flection  of  reality  means  not  that  correlated  conditions  among  components 
produces  higher  reliabilities,  but  rather  that  the  assumption  of  unconditional 
statistical  Independence  of  failures  wes  not  realistic  to  begin  with. 

It  is  possible,  on  the  other  hand,  to  conceive  of  a  situation  in  which 
tne  environmental  profile  is  altered  by  placing  components  in  a  state  of 
physical  dependence.  The  effects  are  reflected  in  correlated  failures  and 
serial  system  reliability  is  actually  increased  over  what  it  would  be  had 
t  <e  components  not  been  placed  in  such  a  state  of  mutual  phys  cal  dependence. 
iv,,  example  is  the  case  where  two  transistors  are  mounted  on  a  common  heat 
vk.  The  assumption  of  statistical  independence  of  components  under  the 
environmental  profile  would  not  be  justified  and  Equation  (2)  would  have  to 
be  nodified.  Thus,  correlations  among  component  conditions  can  be  present 
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with  or  without  the  assumption  of  their  statistical  independence  under  an 
environmental  profile  but  the  interpretation  of  these  correlations  depends 
upon  which  conditions  hold. 

A  second  order  approximation  of  the  difference  between  the  reliability 
estimates  given  by  Equation  (2)  and  the  model  R(t)  ■  (p(t))‘  can  be  obtained 
by  the  use  of  Taylor's  series.  When  the  environment  is  specified  in  terms  of 
a  scalar  x  the  difference  is  calculated  as  follows  for  an  N  component  serial 

system. 

Let  R^Ujx)  -  (R(t,x))W  and  let  w(x)  be  the  probability  density  of  the 
environmental  variable  x.  Equation  (2)  gives  the  system  reliability  as 

RjjU)  -  £  <yt,x»  -  E  [(R(t,x))N]. 


The  Taylors  series  approximation  is 

o 

•  N 

r»  /*.\  _  /n/*.  *  Q 


i  (R(t,x)  )  +  -5-  -S-  (R(tfx)j)N- 

x-U  * X»U 


R‘V,U>  + N.RN’Z(t,yM(N-l)(R'(t,p))2  +  R(t,y)R"(t,y)  J 


where  primes  denote  differentiation  with  respect  to  x. 

The  difference  A  in  the  estimates  is 
2 

A  •  — .N.RN"2(t,y).l(N-l)(R’(tty))2  +  R(t ,y)R"(t,y) ]  (4) 

2 

PROOF  OF  AN  IMPORTANT  INEQUALITY 


The  conclusions  drawn  in  the  previous  section  depend  in  part  upon  the 
validity  of  Equation  (3).  This  inequality  for  the  general  case  follows  from: 


Theorem  2  Let  n  *  (n^»  •••  »  be  a  vector  random  variables  having 
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fj.ni.te  moments  and  having  a  joint  density  function  w(n).  n  is  assumed  to 
define  an  environmental  profile  go  .  Define  N  characteristic  random  variables 
as  follows: 

|  1  if  i-th  component  survives  in  £  for  time  t. 

Xi“  1 

(  0  otherwise 

Let  the  X^'s  be  mutually  independent  under  5  . 

Let  RN(t)  ■  Pr  (N  components  survive  £  for  time  t) . 

Then 

R+k(t)  -  R^t)  R^t)  for  h,  k  positive  integers. 

Proof  By  definition, 

R^t)  -  E(X1  X2  ...  X^)  •  EE(X1  ...  Xj^/n)  -  E{E(X1/n)E(X2/n)...E(X/n)} 

-  E(RN(t,n))  -  E  aN  Where  a  -  R(t,n)  -  E(X  /n) 

n  a 

(0  *  a  ■  1). 

Thus  Rj^  is  the  N-th  raw  moment  of  some  distribution  defined  on  [0,1], 

For  any  positive  integers  h  and  k  and  Independent  random  variables  and 
«2  defined  on  [0,1]  and  having  the  distribution  of  a,  it  follows  that 

(^h  -  -  a2k>  5  0 

so  that 

E(axh+k)  -  2E(a1h  a2k)  +  E(a2h+k)  - 
2E(ah+k)  -  2E(ah)E(ak)  *  0  . 

Thus  L(a  )  *  E(a^)E(ak). 

Rewriting  the  last  inequality  in  terms  of  reliabilities 


V(t)  *  V‘>  Vc> 


iisary  et  al,  proved  a  theorem  analogous  to  Theorem  2  in  terns  of  associated 

random  va'ir.'jl-:-.  n  ,  ...  ,  as  veil  as  proving  another  Inequality  of 
f 'inter  generality. 
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CHAPTER  III 


A  STRESS-STRENGTH  MODEL  OF  RELIABILITY 

The  environmental  variables  model  has  appeared  in  the  literature  in  a 
somewhat  different  form,  Lloyd  and  Lipow  [/»)  os  well  as  others  have  con¬ 
sidered  the  "chain"  or  "weakest  link"  model,  also  called  the  "stress-strength" 
model,  A  generalized  version  of  this  model  for  the  reliability  of  an 
n "Component  serial  system  follows. 

Components  are  assumed  to  function  in  a  random  environment  characterized 
by  a  sequence  of  stress  "pulses",  random  in  number  over  a  time  interval  of 
fixed  length  and  each  following  some  distribution  of  stress  intensity.  Thus 
the  environment  is  defined  in  terns  of  a  discrete  variable, 

Jl(t)  ■  number  of  stress  pulses  in  a  time  interval  of 
length  t 

and  a  distribution  of  stress  Intensity 

F(x)  -  Pr  (a  randomly  selected  stress  pulse  does  not 
exceed  x  in  intensity) 

Each  component  is  assumed  to  possess  a  "strength"  selected  at  random 

from  a  distribution  G(x).  Reliable  operation  of  the n  -component  system  over 

a  time  interval  of  length  t  occurs  whenever  the  minimum  component  strength 

exceeds  the  maximum  stress  intensity  occurring  during  the  period. 

Let  Y  denote  the  maximum  stress  intensity  in  a  time  interval  of  length 
n 

t  given  that  n  stress  pulses  occur. 

Then,  assuming  independence  of  the  pulse  intensities, 

l’(Yn  ■  x)  •  (F(x))°  (n  -  0,  1,  2,  ...  ) 
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Assuming  N(c)  and  Yr  to  be  statistically  independent, 

P(N(t)  »  n,  Y  -  x)  -  P(N(t)  -  n)(F(x))“  (n  -  0,  1,  ...) 
n 

The  marginal  distribution  of  the  maximum  stress  intensity  in  a  time  Interval 
of  length  t,  Y(t) f  is  therefore 

m 

P(Y(t)  *  x)  -  l  P(N(t)  -n)(F(x))n  . 
n»0 

For  example,  assume  that  (N(t);  t  *  0}  is  a  stationary  Poisson  process  of 
Intensity  X. 

Then 

P(Y(t)  *  x)  -  l  e‘Xt(F(x))n 

n«0 


.  e-Xt(l-F(x))  £  (Xt  F(x))n 
n-0  nl 


-  t*F(x)  • 

e 


Thus, 

H(x;t)  •  P(Y(t)  m  x) 


-Xt(l-F(x)) 

e 


(6) 


If  each  component  etrength  follows  the  distribution  G(x)  and  component 
strengths  are  independent,  the  distribution  function  GR(x)  of  the  minimum 
component  strength  X^  is  therefore 

Gn(x)  -  1  -  P(X(n)>  x)  -  1  -  (1  -  G(x))". 

The  n-component  serial  reliability  is  defined  as 

00 

Rn(t)  -  P(X(n)  >  Y(t))  -  /  (1  -  Gn(x))dH(x;t)  (7) 

Thus,  in  the  special  case  given  by  Equation  (6)  and  by  making  the  plausible 
assumption  that  F(0)  ■  0, 

R  (t)  -  e“U  +  /  (1  -  G(x))nd  e-xt(1“F(x))  (7') 

n  Of  • 
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The.  similarity  between  Equation  (7)  and  Equation  (2)  is  obvious. 
Recalling  that  from  Equation  (2)  can  be  interpreted  as  the  n-th  raw 
moment  of  a  dirtrlbution  defined  on  [0,1]  it  is  also  clear  that  R^Ct)  from 
Equation  (7)  can  be  interpreted  in  the  same  manner.  The  model  defined  by 
Equation  (7)  can  be  generalized  to  the  case  of  k  distinct  stress  and  strength 
variables. 
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CHAPTER  IV 


A  MULTIVARIATE  BERNOULLI  MODEL  OF  RELIABILITY 

A  somewhat  different  approach  toward  accounting  for  dependencies  among 
component  conditions  was  taken  by  Patterson  and  Khanna  [5]  in  which  the  basic 
result  is  stated  as 

Theorem  3  Let  X, .....  X  be  n  Bernoulli  random  variables  each  having 
.  1  *  n 

parameter  p  and  possessing  a  joint  probability  function  p(X^,  X2,  ,  X  ) . 

Let  all  pairs  X^,  Xj  have  linear  correlation  P  .  Let  all  conditioned  pairs 

(X.  ,  XA  |x.  -  1,  ,  X,  -  1)  also  have  linear  correlation p  . 

11  2  13  xn 

Then 

n  i-i 

P(X1  -  1 . X^  -  1)  -  it  [1  -  (1-P)  (1-p)]  (8) 

1-1  (0  ip  *  1)  (n  -1,2,3,...) 

If  1,(X1  -  1)  denotes  the  reliability  of  the  i-th  component  then  the  model 
given  by  Equation  (8)  gives  the  n-component  serial  reliability  of  a  system  in 
which  component  conditions  exhibit  a  statistical  dependence.  The  reliability  of 
an  n  component  parallel  system  can  be  obtained  from  (8)  by  the  expression 

1-  P(Xx  -  0,  X2  -  0,  ...  ,  Xn  -  0)  -  IT  [1  -(1-p)  “  p]  (8') 

From  Equation  (8)  it  is  seen  that 

R?  ^  R  ■  R.  (0  *  p  *  1) 
i  n  i 

wherp  R  ■  P(X.  ■  1.  ...  .X  ■  1). 
n  1  n 
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Geometrically  the  probabilities  PCX^,  ...  ,  Xq)  are  mis  assignments  to  the 
vertices  of  an  n-dimenslonal  unit  cube  that  sum  to  unity.  For  the  n  ■  2 
case  Equations  (8)  and  (2)  are  equivalent  if 


P(Xi  -  1)  -  /  R(t,X)uOOdX  <i  -  1,  2). 

3? 

Nothing  is  assumed,  however,  in  the  derivation  of  Equation  (8)  about  the 
cause  mechanism  either  with  respect  to  the  environment  or  possible  physical 
connections  smong  the  components  that  have  the  affect  of  producing  correlations 
between  component  failures.  The  assumption  is,  simply,  that  it  is  possible  to 
express  total  serial  system  reliability  as  a  product  of  conditional  probabilities 
such  that  the  probability  of  survival  of  a  given  component  increases  when  it  is 
known  that  increasing  numbers  of  the  other  n  -  1  components  of  the  system 
survived.  Since  each  conditional  probability  is  of  the  form 


P(X<  -  l|x.  «  1,  X4  -  1,  ...  ,  X*  -  1)  -  1  -  (1-p)  q 
xk+l  11  x2  k 

the  additional  knowledge  that  the  (k+1)  at  component  survived  increases  the 

conditional  probability  that  one  of  the  remaining  n-(lctl)  components  survived 

»k 

by  the  amount  pq(l-p)  . 

Although  the  multivariate  Bernoulli  model  given  by  Equation  (8)  represents 
a  different  approach  to  the  problem  the  sequence 


{r  }  ■  {  it  u-a1’^))  (a  ■  i-o )  ,  oMi  i) 

n  i-i 

(p  ■  l-q) (n  •  1,2,  ...  ) 


represents  the  momenta  of  a  probability  distribution  defined  over  the  closed 
unit  interval.  This  fact  was  demonstrated  in  a  corollary  proved  by  John  Saw 
which  depends  upon  a  theorem  due  to  7.  Hausdorff  and  discussed  by  Feller  [6]. 
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The  statements  of  the  theorem  and  corollary  follow. 


Definition  Let  F  be  a  probability  distribution  concentrated  on  the  interval 
[0,1].  The  k-th  moment  yk  of  F  is 


U.  -  E(Xk)  -  /  Xk  dF. 
K  0 


Ut  E(X  (1  -  X))  -  -  Awk. 
Then  by  induction 


<-A)&k  -  E(Xk(l-X)r). 


Theorem 

A  sequence  of  numbers  y'g,  y^»  •••  represents  the  moments  yk  of  some 
of  some  probability  distribution  F  concentrated  on  [0,1]  if  end  only  if 
(“d)5k  *  0  ,  yQ  ■  1. 

Corollery 

Let  0  <  Uj^  4  u2  <  ...  «  1  be  a  nondecreasing  set  of  positive  values 
bounded  above  by  unity  and  satisfying 


<-«uk  « 


and  define  yA  ■  1  *  y.  ■  *  u  . 


Then  y0,  y^,  ...  are  the  moments  of  some  probability  distribution  concentrated 
on  [0,1]. 

k  i-1 

By  letting  yQ  -  1,  y.  -  d  -  6  q) 


so  that 

i 

H 

a 

o 

ft1"1*, 

p  q  » 

(i 

then 

0  <  \x^  * 

U£  ^  ... 

*  1 

r 

k-1 . 

-  e) 

and 

(-d)u,.  " 

-  6  (1 

The  conditions  of  the  corollary  are  satisfied  and  therefore  the  sequence 


1.  (B  -  1  -  p  -  0) 

For  B  -  0,  RQ  ■  1,  8^  ■  p 
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<k  -  1) 


and 


Than 


and 


P« 


0 


(0  <  k  <  n) 


(q  { 0  *  x  <  1 
lix  -  1 


2.  (e  -  i) 


so  Chat 


and 


\  •  •  kt°. 
p‘")  •  <j>  pk  «  -  p)n‘k 

F  Ck>  •  I  <J)  pk  O  -  P)n  k 
B  lAm 


!  P  ut“»-k  g  ■  du 

o  B([nx] ;  n  -[nx]  ♦  1) 


where  [nx]  la  the  largest  Integer  satisfying  [nx]  <  nx  and  B(r,s)  is  the  Beta 
function  with  arguments  r  and  a. 


Thus 


F(x)  - 


0  l  x  <  p 
i ;  x  -  p 


The  preceding  theory  ahowa  that  the  given  by  Equations  (2) ,  (7) ,  and 
(8)  can  each  be  Interpreted  as  the  n-th  raw  noaenis  of  a  probability  dlstrlbu- 
tion  defined  on  [0,1]. 
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CONCLUSIONS 


The  models  given  by  Equations  (2),  (7),  and  (3)  represent  three,  retlot'. 

of  eli...I  natinp,  the  necessity  of  the  assumption  of  unconditional  statistical 

independence  among  components  operating  jointly  in  a  common  environment.  It 

has  been  shown  that  R  can  be  Interpreted  as  the  n-th  raw  moment  of  a 

n 

probability  distribution  F(x)  distributed  over  [0,1].  It  was  showr  that  the 
assumption  of  the  unconditional  statistical  Independence  of  component  con¬ 
ditions  results  in  an  underestimate  of  serial  system  reliability  and  an  over¬ 
estimate  of  the  reliability  of  a  system  of  components  in  parallel. 

The.  environmental  profile  vector  X  •  (Xi,  . ..,  Xk)  constitutes  a  set  of 
associated  random  variables  whenever  the  reliabilities  p^(X)  and  change 

monot  mically  with  respect  to  X. 

Finally,  the  environmental  profile  and  Bernoulli  models  represent  what 
might  properly  be  called  static  models  in  the  sense  '.hat  none  of  the  random 
varitales  involved  represent  stochastic  processes.  The  stress-strength  node! 
is  d\  nanic  so  far  as  the  representation  of  stress  is  concerned. 


The  author  expresses  his  gratitude  to  John  Saw  for  helpful  criticisms 
and  suggestions  in  the  preparation  of  this  paper. 
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Standard  reliability  prediction  formulas  for  multi-component  systems  make  the- 
assumption  of  statistical  independence  of  the  conditions  of  the  components  after  a 
fixed  period  of  time  in  a  hypothesized  environment.  Although  laboratory  and  field 
experience  shows  that  this  assumption  is  not  always  valid  it  persists  as  a  basis 
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inathematics  of  probability  and  one  soon  discovers  that  the  assumption  of  statistical 
independence  among  components  usually  leads  to  the  simplest  mathematics.  A  second 
reason  why  the  independence  assumption  is  so  predominant  is  that  engineers  and 
analysts  are  not  clear  on  what  alternatives  should  be  pursued  or  even  what  the 
alternatives  might  be  so  there  would  be  little  point  in  making  assumptions  the 
modeling  implications  of  which  are  simply  not  understood..  There  is  another  justifi- 
cantion  for  the  statistical  independence  assumption,  that  being  the  fact  that  such 
ar  assumption  yields  models  that  can  provide  bounds  on  system  reliabilities.  Ana¬ 
lytical  models  which  attempt  to  account  for  environmental  effects  on  component 
failure  rates  show  that  the  statistical  independence  assumption  often  leads  to  gross 
overestimates  or  underestimates  of  system  reliability.  Evans  [1]  argues  that  it  is 
almost  never  the  case  that  this  assumption  is  correct,  but  also  points  out  that 
other  models  which  attempt  to  quantify  or  correct  this  error  have  their  own  limita¬ 
tions.  If,  on  the  other  hand,  it  could  be  shown  that  an  alternative  model  even 
though  an  approximation  is  indeed  a  closer  representation  to  the  true  state  of  the 
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Abstract  (continued) 

s'orld  and  is  at  the  same  time  practical  to  use  in  a  computational  sense  then  it 
should  be  used.  The  question  of  validity  lies  with  the  basic  assumptions  under¬ 
lying  the  model  and  not  with  the  mathematics  itself.  If  a  set  of  assumptions 
alternative  to  the  assumption  of  statistical  independence  are  accepted  as  more 
closely  characterizing  an  environmental  situation  then  the  resulting  mathematical 
predictions  would  be  tentatively  preferred,  subject  to  experimental  verification. 

This  paper  reviews  some  previous  explorations  of  the  question  of  statistical 
dependence  of  operating  components  and  the  effects  on  system  reliability  and 
presents  some  results  not  previously  worked  out.  In  particular  there  is  shown 
co  exist  a  connection  among  three  models  that  do  not  require  the  assumption  of 
ctatistical  independence  of  components.  Under  certain  conditions  these  models 
ran  generate  the  same  reliability  prediction  even  though  they  appear  to  be 
drived  under  different  assumptions. 


